A framework for real-time monitoring, analysis and adaptive sampling of viral amplicon nanopore sequencing.

The ongoing SARS-CoV-2 pandemic demonstrates the utility of real-time sequence analysis in monitoring and surveillance of pathogens. However, cost-effective sequencing requires that samples be PCR amplified and multiplexed via barcoding onto a single flow cell, resulting in challenges with maximising and balancing coverage for each sample. To address this, we developed a real-time analysis pipeline to maximise flow cell performance and optimise sequencing time and costs for any amplicon based sequencing. We extended our nanopore analysis platform MinoTour to incorporate ARTIC network bioinformatics analysis pipelines. MinoTour predicts which samples will reach sufficient coverage for downstream analysis and runs the ARTIC networks Medaka pipeline once sufficient coverage has been reached. We show that stopping a viral sequencing run earlier, at the point that sufficient data has become available, has no negative effect on subsequent down-stream analysis. A separate tool, SwordFish, is used to automate adaptive sampling on Nanopore sequencers during the sequencing run. This enables normalisation of coverage both within (amplicons) and between samples (barcodes) on barcoded sequencing runs. We show that this process enriches under-represented samples and amplicons in a library as well as reducing the time taken to obtain complete genomes without affecting the consensus sequence.

Spike Mutation Profiles Associated With SARS-CoV-2 Breakthrough Infections in Delta Emerging and Predominant Time Periods in British Columbia, Canada

Background COVID-19 vaccination is a key public health measure in the pandemic response. The rapid evolution of SARS-CoV-2 variants introduce new groups of spike protein mutations. These new mutations are thought to aid in the evasion of vaccine-induced immunity and render vaccines less effective. However, not all spike mutations contribute equally to vaccine escape. Previous studies associate mutations with vaccine breakthrough infections (BTI), but information at the population level remains scarce. We aimed to identify spike mutations associated with SARS-CoV-2 vaccine BTI in a community setting during the emergence and predominance of the Delta-variant. Methods This case-control study used both genomic, and epidemiological data from a provincial COVID-19 surveillance program. Analyses were stratified into two periods approximating the emergence and predominance of the Delta-variant, and restricted to primary SARS-CoV-2 infections from either unvaccinated individuals, or those infected ≥14 days after their second vaccination dose in a community setting. Each sample's spike mutations were concatenated into a unique spike mutation profile (SMP). Penalized logistic regression was used to identify spike mutations and SMPs associated with SARS-CoV-2 vaccine BTI in both time periods. Results and Discussion This study reports population level relative risk estimates, between 2 and 4-folds, of spike mutation profiles associated with BTI during the emergence and predominance of the Delta-variant, which comprised 19,624 and 17,331 observations, respectively. The identified mutations cover multiple spike domains including the N-terminal domain (NTD), receptor binding domain (RBD), S1/S2 cleavage region, fusion peptide and heptad regions. Mutations in these different regions imply various mechanisms contribute to vaccine escape. Our profiling method identifies naturally occurring spike mutations associated with BTI, and can be applied to emerging SARS-CoV-2 variants with novel groups of spike mutations.

Two-Dose Severe Acute Respiratory Syndrome Coronavirus 2 Vaccine Effectiveness With Mixed Schedules and Extended Dosing Intervals: Test-Negative Design Studies From British Columbia and Quebec, Canada.

BACKGROUND: The Canadian coronavirus disease 2019 (COVID-19) immunization strategy deferred second doses and allowed mixed schedules. We compared 2-dose vaccine effectiveness (VE) by vaccine type (mRNA and/or ChAdOx1), interval between doses, and time since second dose in 2 of Canada's larger provinces. METHODS: Two-dose VE against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection or hospitalization among adults ≥18 years, including due to Alpha, Gamma, and Delta variants of concern (VOCs), was assessed ≥14 days postvaccination by test-negative design studies separately conducted in British Columbia and Quebec, Canada, between 30 May and 27 November (epi-weeks 22-47) 2021. RESULTS: In both provinces, all homologous or heterologous mRNA and/or ChAdOx1 2-dose schedules were associated with ≥90% reduction in SARS-CoV-2 hospitalization risk for ≥7 months. With slight decline from a peak of >90%, VE against infection was ≥80% for ≥6 months following homologous mRNA vaccination, lower by ∼10% when both doses were ChAdOx1 but comparably high following heterologous ChAdOx1 + mRNA receipt. Findings were similar by age group, sex, and VOC. VE was significantly higher with longer 7-8-week versus manufacturer-specified 3-4-week intervals between mRNA doses. CONCLUSIONS: Two doses of any mRNA and/or ChAdOx1 combination gave substantial and sustained protection against SARS-CoV-2 hospitalization, spanning Delta-dominant circulation. ChAdOx1 VE against infection was improved by heterologous mRNA series completion. A 7-8-week interval between first and second doses improved mRNA VE and may be the optimal schedule outside periods of intense epidemic surge. Findings support interchangeability and extended intervals between SARS-CoV-2 vaccine doses, with potential global implications for low-coverage areas and, going forward, for children.

Comparative Single-Dose mRNA and ChAdOx1 Vaccine Effectiveness Against Severe Acute Respiratory Syndrome Coronavirus 2, Including Variants of Concern: Test-Negative Design, British Columbia, Canada.

BACKGROUND: In British Columbia, Canada, most adults 50-69 years old became eligible for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccine in April 2021, with chimpanzee adenoviral vectored vaccine (ChAdOx1) restricted to ≥55-year-olds and second doses deferred ≥6 weeks to optimize single-dose coverage. METHODS: Among adults 50-69 years old, single-dose messenger RNA (mRNA) and ChAdOx1 vaccine effectiveness (VE) against SARS-CoV-2 infection and hospitalization, including variant-specific, was assessed by test-negative design between 4 April and 2 October 2021. RESULTS: Single-dose VE included 11 861 cases and 99 544 controls. Median of postvaccination follow-up was 32 days (interquartile range, 15-52 days). Alpha, Gamma, and Delta variants comprised 23%, 18%, and 56%, respectively, of genetically characterized viruses. At 21-55 days postvaccination, single-dose mRNA and ChAdOx1 VE (95% confidence interval [CI]) was 74% (71%-76%) and 59% (53%-65%) against any infection and 86% (80%-90%) and 94% (85%-97%) against hospitalization, respectively. VE (95% CI) was similar against Alpha and Gamma infections for mRNA (80% [76%-84%] and 80% [75%-84%], respectively) and ChAdOx1 (69% [60%-76%] and 66% [56%-73%], respectively). mRNA VE was lower at 63% (95% CI, 56%-69%) against Delta but 85% (95% CI, 71%-92%) against Delta-associated hospitalization (nonestimable for ChAdOx1). CONCLUSIONS: A single mRNA or ChAdOx1 vaccine dose gave important protection against SARS-CoV-2, including early variants of concern. ChAdOx1 VE was lower against infection, but 1 dose of either vaccine reduced the hospitalization risk by >85% to at least 8 weeks postvaccination. Findings inform program options, including longer dosing intervals.

Single-dose mRNA Vaccine Effectiveness Against Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), Including Alpha and Gamma Variants: A Test-negative Design in Adults 70 Years and Older in British Columbia, Canada.

BACKGROUND: Randomized-controlled trials of messenger RNA (mRNA) vaccine protection against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) included relatively few elderly participants. We assess single-dose mRNA vaccine effectiveness (VE) in adults ≥ 70 years old in British Columbia, Canada, where second doses were deferred by up to 16 weeks and where a spring 2021 wave uniquely included codominant circulation of Alpha (B.1.1.7) and Gamma (P.1) variants of concern (VOC). METHODS: Analyses included community-dwelling adults ≥ 70 years old with specimen collection between 4 April (epidemiological week 14) and 1 May (week 17) 2021. Adjusted VE was estimated by test-negative design. Cases were reverse-transcription polymerase chain reaction (RT-PCR) test-positive for SARS-CoV-2, and controls were test-negative. Vaccine status was defined by receipt of a single-dose ≥ 21 days before specimen collection, but a range of intervals was assessed. Variant-specific VE was estimated against viruses genetically characterized as Alpha, Gamma or non-VOC lineages. RESULTS: VE analyses included 16 993 specimens: 1226 (7%) test-positive cases and 15 767 test-negative controls. Of 1131 (92%) genetically characterized viruses, 509 (45%), 314 (28%), and 276 (24%) were Alpha, Gamma, and non-VOC lineages, respectively. At 0-13 days postvaccination, VE was negligible at 14% (95% confidence interval [CI], 0-26) but increased from 43% (95% CI, 30-53) at 14-20 days to 75% (95% CI, 63-83) at 35-41 days postvaccination. VE at ≥ 21 days postvaccination was 65% (95% CI, 58-71) overall: 72% (95% CI, 58-81), 67% (95% CI, 57-75), and 61% (95% CI, 45-72) for non-VOC, Alpha, and Gamma variants, respectively. CONCLUSIONS: A single dose of mRNA vaccine reduced the risk of SARS-CoV-2 by about two-thirds in adults ≥ 70 years old, with protection only minimally reduced against Alpha and Gamma variants.

Mutations in emerging variant of concern lineages disrupt genomic sequencing of SARS-CoV-2 clinical specimens.

Mutations in emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages can interfere with laboratory methods used to generate viral genome sequences for public health surveillance. We identified 20 mutations that are widespread in variant of concern lineages and affect widely used sequencing protocols by the ARTIC network and Freed et al. Three of these mutations disrupted sequencing of P.1 lineage specimens during a recent outbreak in British Columbia, Canada. We provide laboratory validation of protocol modifications that restored sequencing performance. The study findings indicate that genomic sequencing protocols require immediate updating to address emerging mutations. This work also suggests that routine monitoring and protocol updates will be necessary as SARS-CoV-2 continues to evolve. The bioinformatic and laboratory approaches used here provide guidance for this kind of assay maintenance.

Assessing Multiplex Tiling PCR Sequencing Approaches for Detecting Genomic Variants of SARS-CoV-2 in Municipal Wastewater.

Wastewater-based genomic surveillance of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus shows promise to complement genomic epidemiology efforts. Multiplex tiling PCR is a desirable approach for targeted genome sequencing of SARS-CoV-2 in wastewater due to its low cost and rapid turnaround time. However, it is not clear how different multiplex tiling PCR primer schemes or wastewater sample matrices impact the resulting SARS-CoV-2 genome coverage. The objective of this work was to assess the performance of three different multiplex primer schemes, consisting of 150-bp, 400-bp, and 1,200-bp amplicons, as well as two wastewater sample matrices, influent wastewater and primary sludge, for targeted genome sequencing of SARS-CoV-2. Wastewater samples were collected weekly from five municipal wastewater treatment plants (WWTPs) in the Metro Vancouver region of British Columbia, Canada during a period of increased coronavirus disease 19 (COVID-19) case counts from February to April 2021. RNA extracted from clarified influent wastewater provided significantly higher genome coverage (breadth and median depth) than primary sludge samples across all primer schemes. Shorter amplicons appeared to be more resilient to sample RNA degradation but were hindered by greater primer pool complexity in the 150-bp scheme. The identified optimal primer scheme (400 bp) and sample matrix (influent) were capable of detecting the emergence of mutations associated with genomic variants of concern, for which the daily wastewater load significantly correlated with clinical case counts. Taken together, these results provide guidance on best practices for implementing wastewater-based genomic surveillance and demonstrate its ability to inform epidemiology efforts by detecting genomic variants of concern circulating within a geographic region. IMPORTANCE Monitoring the genomic characteristics of the SARS-CoV-2 virus circulating in a population can shed important insights into epidemiological aspects of the COVID-19 outbreak. Sequencing every clinical patient sample in a highly populous area is a difficult feat, and thus sequencing SARS-CoV-2 RNA in municipal wastewater offers great promise to augment genomic surveillance by characterizing a pooled population sample matrix, particularly during an escalating outbreak. Here, we assess different approaches and sample matrices for rapid targeted genome sequencing of SARS-CoV-2 in municipal wastewater. We demonstrate that the optimal approach is capable of detecting the emergence of SARS-CoV-2 genomic variants of concern, with strong correlations to clinical case data in the province of British Columbia. These results provide guidance on best practices on, as well as further support for, the application of wastewater genomic surveillance as a tool to augment current genomic epidemiology efforts.

Rapid Increase in SARS-CoV-2 P.1 Lineage Leading to Codominance with B.1.1.7 Lineage, British Columbia, Canada, January-April 2021.

Several severe acute respiratory syndrome coronavirus 2 variants of concern (VOCs) emerged in late 2020; lineage B.1.1.7 initially dominated globally. However, lineages B.1.351 and P.1 represent potentially greater risk for transmission and immune escape. In British Columbia, Canada, B.1.1.7 and B.1.351 were first identified in December 2020 and P.1 in February 2021. We combined quantitative PCR and whole-genome sequencing to assess relative contribution of VOCs in nearly 67,000 infections during the first 16 weeks of 2021 in British Columbia. B.1.1.7 accounted for <10% of screened or sequenced specimens early on, increasing to >50% by week 8. P.1 accounted for <10% until week 10, increased rapidly to peak at week 12, and by week 13 codominated within 10% of rates of B.1.1.7. B.1.351 was a minority throughout. This rapid expansion of P.1 but suppression of B.1.351 expands our understanding of population-level VOC patterns and might provide clues to fitness determinants for emerging VOCs.